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Abstract 

We consider the genealogy tree for a critical branching process con- 
ditioned on non-extinction. We enumerate vertices in each generation 
of the tree so that for each two generations one can define a mono- 
tone map describing the ancestor-descendant relation between their 
vertices. We show that under appropriate rescaling this family of 
monotone maps converges in distribution in a special topology to a 
limiting flow of discontinuous monotone maps which can be seen as a 
continuum tree. This flow is a solution of an SPDE with respect to a 
Brownian sheet. 
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1 Introduction 

In this paper we present a point of view at large random trees. In [Bak] 
we considered Boltzmann-Gibbs distributions on rooted plane trees with 
bounded branching and proved that as the order of the tree grows to in- 
finity, these trees obey a certain thermodynamic limit theorem with a limit 
given by an infinite Markov random tree that can be seen as the genealogy 
tree for a specially chosen critical branching processes conditioned on non- 
extinction, see [AP98j,[Kes86],[Ken75]. Appropriately rescaled generation 
sizes in the limiting tree admit an approximation by a diffusion process. 
However, this view ignores all the interesting details concerning the com- 
plicated way different generations of the random tree are connected to each 
other. 

The goal of this paper is to extend the diffusion approximation result 
and investigate the fine structure of the infinite random tree. The idea is 
to encode the tree via a stochastic flow of monotone maps. We show that 
under the same rescaling, the flow of monotone maps associated to the tree 
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converges in distribution to a monotone flow that can be viewed as a solution 
of an evolutionary SPDE with respect to a Brownian sheet. The possibility 
of such a result was hinted at in [Bak] , and in this text we provide a rigorous 
treatment of the SPDE, some interesting properties of its solutions, and the 
convergence theorem in an appropriate metric space. The solution of SPDE 
that we construct can be called a continuum random tree, and it appears 
to be a new type of continuum random trees, unknown in the literature, 
although some other stochastic flows similar to ours have appeared in the 
literature (see, e.g., }LG99| ) and there is a natural connection of our results 
to superprocesses that describe the evolution of mass generated by continual 
branching particle systems conditioned on nonextinction, see [EP90],[Eva93 . 
We describe this connection and the advantages of our approach in Section [9j 
In short, the results of this paper along with those of [Bak] show what a 
typical large tree looks like if drawn on the plane. 

Acknowledgment. The author wishes to thank Tom Kurtz, Carl 
Mueller and Andrey Pilipenko for their comments on connection of this 
work to the existing literature on superprocesses. He is grateful to NSF for 
partial support through CAREER grant DMS-0742424. 

2 Gibbs distributions on plane trees and thermo- 
dynamic limit 

In this section we describe the background in detail. Recall that plane (or, 
ordered) trees are rooted trees such that subtrees at any vertex are linearly 
ordered. 

We fix D G N and introduce = Tn(D), the set of all plane trees 
on iV vertices such that the branching number (i.e. the number of children) 
of each vertex does not exceed D. 

We assume that every admissible branching number i € {0, . . . , D} is 
assigned an energy E^ £ M, and the energy of the tree T € Yn is defined via 

D 

E(T)= Y, E deg{v) =Y,Xi{T)E u 
ueV(T) j=o 

where V(T) denotes the set of vertices of the tree T, deg(v) denotes the 
branching number of vertex v, and Xi(T) is the number of vertices of degree i 
in T. This energy function defines a local interaction between the vertices 
of the tree since the energy contribution from a single vertex depends only 
on its nearest neighborhood in the tree. 
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We fix an inverse temperature parameter f5 > and define a probability 
measure /Ujv on Tjv by 



In [Bakj . the limiting behavior of measures pn as N — > oo was studied. 
To recall the results of [Bakj , we need more notation and terminology. 

For each vertex v of a tree T € Tjv its height h(v) is defined as the 
distance to the root of T, i.e. the length of the shortest path connecting v 
to the root along the edges of T. The height of a finite tree is the maximum 
height of its vertices. 

Each rooted plane tree can be uniquely encoded as a sequence of gen- 
erations. By a generation we mean a monotone (nondecreasing) map G : 
{1, . . . , k} — > N, or, equivalently, the set of pairs {(i, G(i)) : i = 1, . . . , k} 
such that if i\ < i 2 then G{i\) < Gfo)- We denote \G\ = k the number of 
vertices in the generation, and for any % = 1, . . . , k, G(i) denotes i's parent 
number in the previous generation. 

For two generations G and G' we write G < G 1 and say that G' is a 
continuation of G if G'(\G'\) < \G\. Each tree of height n can be viewed as 
a sequence of generations 



where 1 < G\ means G\ i\G\ \ ) = 1 (the 0-th generation consists of a unique 
vertex, the root) and G n < means that the generation n + 1 is empty. 
Infinite sequences 1 < G\ < G2 < • • ■ naturally encode infinite trees. 

For any plane tree T and any n € N, 7r n T denotes the neighborhood of 
the root of radius n, i.e. the subtree of T spanned by all vertices with height 
not exceeding n. 

For any n and sufficiently large N, the map TT n pushes the measure /in 
on Tn forward to the measure Hn^ 1 on S n , the set of all trees with height n. 
We introduce p and C as a unique solution of 




Z N = J2 e~ mT) - 



TeT]v 



1 < Gi < G 2 < ■ ■ ■ < G n < 




and define 



p i = Ce- pEi p i , i = 0, 1,...,D. 
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Then the vector (pi)f=o defines a probability distribution on {0, 1, . . . ,D} 
with mean 1. This vector plays the role of a minimizer for the free en- 
ergy (or large deviation rate function) associated to this model, sec [BH08J 
and [BH09j . 

Some of the results of |Bak] are summarized in the next theorem. 

Theorem 1 1. There is a unique measure P on infinite rooted plane 
trees with branching bounded by D such that for any n £ N 

^7V7T n -> Pn n . 

2. Measure P defines a Markov chain on generations (G n ) n >Q. The tran- 
sition probabilities are given by 

P{G n+1 =g'\G n = g}={^--- P ^ 9<9 '> (1) 

[0, otherwise, 

where i\~, k = 1, . . . , \g\ denotes the number of vertices in generation 
g' that are children of k-th vertex in generation g. 

Remark 1 Our notation differs from the notation used in jBakJ. It is an 
easy exercise to check that the r.h.s. of (H|) equals the expression given 
in JBakJ. 

Remark 2 The infinite Markov tree from Theorem [7] can be obtained as a 
genealogy tree of a critical branching process conditioned on nonextinction, 
see \AP98i . \Kes86j . and \Kenl3j . 

Let us also recall a functional limit theorem for X n = |G n |, the process 
of the Markov infinite tree's generation sizes. We introduce moments of the 
distribution p: 

D 

B n = Y,i n Pi, rien, (2) 
i=0 

and its variance 

H = B 2 - B\ = B 2 - 1. 
Let us now fix a positive time T and define 

Z n it) = —(X[ nt ] + {nt}{X [nt+1] - X [nt] )), t G [0,T], 

the rescaled process of linear interpolation between values of Xk given at 
whole times k = 0, 1, 
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Theorem 2 The distribution of Z n converges weakly in uniform topology 
on C([0, T]) to the distribution of a diffusion process with generator 

Lf(x) = (f'(x) + ^xf"(x))l {xm , (3) 
emitted from at time 0. 

Notice that the limiting process can be viewed as a weak solution of an SDE 

dX{t) = dt + ^/X{t)dW{t) 

on R + (zero is entrance and non-exit singular boundary point for K+). 

This result captures the asymptotics of a rather rough characteristic of 
the tree, the generation size. It ignores all the interesting details of the 
way the generations are connected to each other. In this paper, we are 
interested in taking these details into account, thus refining the diffusion 
approximation. 

Prom now on we consider the random genealogy trees defined above via 
the probability vector (pj)^_ - Although D will always be assumed to be 
finite, the results also hold true if one replaces the boundedness of D with 
certain moment restrictions on the distribution. 



3 Limit theorems for finite partitions 

Let us take a positive number z, and for each n £ N consider the Markov 
process on generations (described in the previous section) originating at 
time with a population of [nz] vertices. Let us take m € N and choose 
some partition = xq < x\ < . . . < x m _i < x m = z of [0,z]. Let us de- 
note by U(n, k,j) the size of the progeny of first [nx^] vertices in the initial 
population after j steps. Notice that V(n, k, j) = U(n,k,j) — U(n,k — 1, j) 
describes the evolution of the progeny of disjoint subpopulations of the orig- 
inal population. Denoting the rescaled interpolation by 

Vn,k{t) = —{V(n,k,[nt])+{nt}(V(n,k,[nt] + l)-V(n,k,[nt]))), t G [0, T], 
/in 

and following the same lines as in |Bak| . we obtain the following result: 

Theorem 3 Vector-valued process (V n< k)^ =1 converges (as n — > oo) in dis- 
tribution in the uniform topology to a diffusion process (Vfc)P on M.™ with 
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initial data 14(0) = v k = x k — x k -\ and generator L, for all C 2 -functions 
with compact support given by 

m 1 m 

Lf(y) = Y,^ 1 {y k >o}dkf(y) + ^J2y^ 1 {y k >o}dl k f(y), ^K™ (4) 
fc=l mi 1 fe=i 

where \y\\ =yi + ... + y m . 

Equivalently, the limiting process can be represented as a coordinatewise 
nonnegative solution of a system of SDEs: 

dV ^) = W^T dt + VW)dW k (t), (5) 
V k (0)=v k , k = l,...,m, (6) 

where {W k )' k n =1 are independent Wiener processes. 

It is clear that this system has a unique strong solution until the first 
time at which V k = for some k. It is also easy to see (essentially from 
the Feller classification of singular points for 1-dimensional SDEs) that for 
each k = 1, . . . , m, the coordinate hyperplane {x k = 0} is an absorbing set. 
Therefore, as soon as the system reaches one of the coordinate planes, we 
have to solve a system of m — 1 equations from that point on. Iterating this 
process we obtain a unique strong solution in M™. 

In the same way, a sequence (C^n,fc)fcLi defined by 

U n , k (t) = —(U(n,k,[nt}) + {nt}(U(n,k,[nt]+l)-U(n,k,[nt}))), t € [0,T], 
///> 

satisfies a functional limit theorem. 

Theorem 4 The vector-valued process (U ritk ) k n =1 converges in distribution 
in uniform topology to a limiting process (U n )' l ^_ 1 that can be represented as 

U k {t) = V 1 (t) + ... + V k (t), k = l,...,m, 

where V is the limiting process from Theorem 

Equivalently, the limiting process can be represented as a (nondecreasing 
in k) solution of a system of SDEs: 



Uk(t) 
U m (t) 

£4(0) = x k , k = l,...,m, (8) 



dU k (t) = jP^dt + VUk(t)-U k -i(t)dW 3 (t), (7) 
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where (Wj)'JL 1 are independent Wiener processes. 

Although the theorems above provide interesting information about the 
behavior of the progeny of finitely many subpopulations, we do not give de- 
tailed proofs of these results, since we are actually aiming at more advanced 
ones (nevertheless, see Section fe.31 for a proof of a similar statement). 

It is important to notice that the system I©,© (or, equivalently, ©,([6])) 
not only determines the evolution of process U (respectively V) in time, but 
also has a rich "spatial" structure which will be explored in later sections. 

4 SPDE approximation: first steps 

In this section we give an informal view at the possibility of constructing a 
limiting continuum random tree represented as a solution to a certain SPDE. 

In the approach described in the last section, for any partition (xi-)V? =0 
we were able to find a probability space and a diffusion process defined on it 
that serves as a weak limit for the evolution of subpopulations in a discrete 
tree. Note that for any t > 0, the pre-limit processes V n ,U n as well as 
their limits U, V constructed above give rise to a random monotone map 
Xk t->- Uk(t) defined on (a;fc)fcLo an d the aforementioned probability space. 

A drawback of this approach is that we need a separate probability space 
for each partition whereas it would be natural to have just one probability 
space and a random flow of monotone maps x t— > U(x,t) (defined on that 
space) that serves all partitions simultaneously. That is, we would like 
to require that for any m and any partition (xfc)fclo of [Q,z], (U(xk,t) — 
U(xk-i,t))™ =1 is a Markov process with generator L defined in (jlj). 

To that end we utilize the spatial structure that has been mentioned in 
the end of Section [3] and introduce a stochastic equation with respect to 
a Brownian sheet W defined on Borel subsets of = {(x,t) : x,t > 0} 
and viewed as an orthogonal martingale measure (see |Wal86l Chapter 2] 
for a definition of Brownian sheet and stochastic integration with respect to 
orthogonal martingale measures): 



dU(x,t) 




dt + W{[0,U(x,t)} x dt) 



U(x,t ) 



x, x £ [0, z] 



or, in integral form 



U(x,t) = x + 



L 



o U(z,s) 



ds + 



L 




l[o,U(x,t)](y)W(dy x dt), (9) 
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The choice of this SPDE is a natural one, since the drift term and the 
quadratic covariation structure of the martingale part coincide with those 
of ©. 

Our discussion of this SPDE will concentrate around the following ques- 
tions: a natural definition of a solution, its existence and uniqueness, regu- 
larity properties, and the relation to the finite-dimensional system (J?]),©. 

First, let us take any z > and x = z in ([9]). The equation rewrites then 
as t 

U(z,t) = z + t+ f f l[o,u(z,t)){vW(dy x dt). (10) 

JO JR 

One can obtain a unique strong solution of this equation using the stan- 
dard Picard iteration scheme and mimicking the proof of existence and 
uniqueness for an ordinary SDE. The only potential difficulty one may en- 
counter is the behavior of the solution near 0, but the solution is a Markov 
process with generator L given in ([3]) so that the solution always stays pos- 
itive with probability 1. 

Next, we can take any z > and a finite set I of admissible initial points, 

i.e., 

{0,z} C/C [0,*]. 

Mimicking the construction and the proof of the uniqueness of a positive 
solution of (J!]),© we see that ([9]) has a unique strong nondecreasing in 
x € I solution (U(x,t)) x£ ij>o- Thus obtained U satisfies a finite system of 
stochastic equations w.r.t. the Brownian sheet with probability one. It is 
clear that if I' D I is a broader finite set of admissible initial conditions, 
then, with probability 1, for any t > 0, the map x h-> U(x,t),x £ I' is a 
monotone extension of the map x h4 U(x,t),x € I. 

In principle, there might be a problem with defining this map with prob- 
ability 1 for an uncountable set x £ [0, z]. However, we can use the mono- 
tonicity to tackle this difficulty. Let us take a countable set J of admissible 
initial conditions such that J is dense in [0,z]. Then, taking a sequence of 
finite sets increasing to J, and iteratively extending the map x U(x, t) on 
each set of the sequence as described above, we are able to define this map 
on J with probability 1, for all t > 0. Moreover, this map is a. s. -monotone 
for all t. Therefore, this map can be extended uniquely by monotonicity to 
all points of [0, z] except for at most countable set of discontinuities. Notice 
that for two different dense sets J, J' of admissible initial conditions, with 
probability 1, the maps U defined on J and J' are monotone continuations 
of each other. Therefore the monotone extensions of these two maps agree 
at all points of [0, z\ except for a countable set of upward jumps. 
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It is noteworthy that discontinuities in the form of upward jumps, or 
shocks, do occur with probability 1 as we shall see in Section [7J To deal 
with them we are going to ignore the value of U(x, t) at the jump, and work 
with equivalence classes of monotone functions coinciding at every continuity 
point. We proceed to introduce these equivalence classes in the next section. 

5 Monotone graphs and flows 

The goal of this section is to introduce a space that will serve as a natural 
existence and uniqueness class of solutions of SPDE ^ . 

Consider all points z > and nonnegative nondecreasing functions / 
defined on (— oo, z] such that f(x) = for all x < 0. Each of these functions 
has at most countably many discontinuities. We say that two such functions 
fx : (— oo, Zx] — > M + , /a : (—00,2:2] — > M + are equivalent if zx = Z2, fi{zx) = 
72(2:2), and for each continuity point x of fx, fx(x) = ftix)- Although the 
roles of fx and f2 seem to be different in this definition, it is easily seen to 
define a true equivalence relation. The set of all classes of equivalence will 
be denoted by M. We would like to endow M with a metric structure, and 
to that end we develop a couple of points of view. 

Sometimes, it is convenient to identify each element of M with its unique 
right-continuous representative. Sometimes, it is also convenient to work 
with graphs. The graph of a monotone function / defined on (—00, z] is the 
set Gf = {(x,/(x)), x < z}. For each discontinuity point x of / one may 
consider the line segment f(x) connecting points [x,f(x—)] and [x,f(x+)]. 
The continuous version of Gf is the union of Gf and all segments f(x). It is 
often convenient to identify an element of M with a continuous version of its 
graph restricted to Rj_, and we shall do so from now on calling the elements 
of M monotone graphs. Yet another way to look at monotone graphs is 
to think of them as monotone multivalued maps so that the image of each 
point is either a point f(x) or a segment f(x). 

The Hausdorff distance between Tx € M and I^ € M is defined via 



Lemma 1 (M, p) is a Polish (complete and separable) metric space. 

Proof: The separability follows since one can easily approximate any 
monotone graph by broken lines with finitely many rational vertices. 

The metric space of compact subsets of the plane is complete with respect 
to Hausdorff metric. Therefore, to establish the completeness of (M, p) we 
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need to prove that a limit of a convergent sequence of monotone graphs is 
necessarily a monotone graph. 

The limit set is obviously a connected one, with at most one point on each 
of the lines {(xo,x±) : u$xq + a\X\ = s} for any a$,a.\ > and any s£t. 
Therefore it is a curve {70 (s), 7i(s)} parametrized by s = 70 (s) + 7i(s), and 
the monotonicity follows easily. □ 

We can derive another view at the space M from the proof above. We can 
look at each monotone graph in M in a rotated coordinate frame: then we 
shall see that the set {(7o(s) +7i(s), 71 (s) — 7o( s ))} is a graph of 1-Lipschitz 
function. 

We can introduce another distance based on this viewpoint. Take any 
F±, T 2 G M and consider these curves in M. 2 ^ in the rotated coordinate system 
as 1-Lipschits functions g± and g 2 over [0, si] and [0, S2] respectively. Define 

= s\ A S2, and 

p'{T 1 ,r 2 ) = \s 1 - s 2 \+ sup - ff2(s) |- 

0<s<s» 

The proof of the following lemma an easy exercise. 

Lemma 2 Distances p and p' are equivalent on M. 

The following is a useful criterion of convergence in (M, p). 

Lemma 3 A sequence of monotone graphs (r n ) ne ^ with right- continuous 
representatives f n : [0,z n ] — > M + ,n G N, converges in p to T with right- 
continuous representative f : [0, z\ — > R + iff 

1. z n — >• z, as n — > 00; 

2. fn(z n ) -> f{z), as 00; 

3. At each continuity point x < z of f, f n (x) — > f{x), as n — > 00. 

Proof: Suppose conditions 1-3 are satisfied. If z = 0, then the convergence 
is obvious. If z > 0, take any £ > and find continuity points x\ < x 2 < 
. . . < x m so that x\ < e, x 2 — x\ <£,..., z — x m < e. Then take no € N 
such that for all n > uq, 

(i) for all k = 1, . . . , m, Xk < z n ; 

(ii) for all k = 1,. . . ,m, \f„(x k ) - f(x k )\ < s; 

(iii) \z - z n \ < e; 
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(iv) \f n (z) -f{z n )\ < e. 

It is easy to see that due to the monotonicity, for these values of n we have 
p(r n ,r) < 2e, and thus the convergence holds. 

Suppose now that the convergence holds. Conditions 1 and 2 of the 
lemma follow immediately by monotonicity. If condition 3 is violated, then 
there is 5 > and a subsequence (n 1 ) such that \f n '( x ) — f(x)\ > 5. We can 
choose £ > so that \y — x\ < e implies \f(x) — f(y)\ < 5/2. Therefore, for 
each n', the point (x, f n >(x)) is at least at distance e A (5/2) from any point 
on r which contradicts the assumption and completes the proof. □ 

The criterion of convergence provided by the above lemma allows to con- 
clude immediately that any bounded set of monotone graphs is precompact. 
Let us introduce Zj(T) = sup{xj : (xo,xi) E T}, j = 0, 1. 

Lemma 4 A set H G M is precompact in M iff the set {(z (T), ^i(r)) : T £ 
H} is bounded in M. 2 . 

This lemma can also be derived from Lemma [2] since the Lipschitz con- 
stant for monotone graphs in the rotated coordinate system is bounded by 1. 

Using Lemma El we can also prove the following statement, which we 
will not use but give here for completeness: 

Lemma 5 For any z > 0, the set M(z) of all monotone graphs T € M 
with zq(T) = z is a closed set. Convergence of monotone graphs in M.(z) 
in distance p is equivalent to essential convergence of their monotone right- 
continuous representatives which is in turn equivalent to their convergence 
in Skorokhod topology on D([0,z]). 

For two monotone graphs T\ and T2 with zi(Ti) = -zo(L2)> we define their 
composition T2 ° Ti as the set of all pairs (xo,xi) such that (xq,X2) € Ti 
and (x2,xi) G T2 for some X2- 

Let T > and A T = {(t ,ti) : < t < h < T}. We say that 
(;po, 1 )( t0) t 1 ) e A r is a (continuous) monotone flow on [0,T] if the following 
properties are satisfied: 

1. For each (to,ti) € At, T to,tl is a monotone graph. 

2. The monotone graph r* 0,il depends on (£0, ii) continuously in p. 

3. For each t € [0, T], T t,i: is the identity map on [0, Z(t)] for some Z(t). 
The function Z is called the profile of T. 
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4. For any (t ,h) G A T , zo^* '* 1 ) = Z(t ), z^T***) = Z(h), where Z 
is the profile of T. 

5. If (i ,ii) G A T and (ti,t 2 ) G A T , then r* '* 2 = r* 1 '* 2 o r* '* 1 . 

It is easy to check that the space M[0, T] of all monotone flows on [0, T] 
is a Polish space if equipped with uniform Hausdorff distance: 

pt(Ti,t 2 )= sup pcrf^if'* 1 ). (ii) 

(t ,*i)eA T 

Property O (consistency) implies that Property [2] (continuity) has to be 
checked only for to = H- 

Our next goal is to introduce trajectories of individual points in the 
monotone flow. Suppose T G At, and Z is the profile of T. Let U : 
M + x Ax — > K + satisfy the following properties: 

1. For any (to,t{) G Ay, the function U(x,to,ti) is monotone in i £ 
[0,Z(t )]. 

2. For any (to,ti) € A T , if x € [0,Z(t )], then (ti,C/"(ar,t ,ti)) G T* '* 1 

3. For all x,to, ^(sc^Oj^i) is continuous in ti. 

Then f7 and T are said to be compatible with each other, and U is said to 
be a trajectory representation of V. Clearly, the monotonicity implies that, 
given Z, there is at most one monotone flow on [0, T] compatible with U. 
Moreover, it is sufficient to know a trajectory representation U(x,to,t\) for 
a dense set of points x,to,ti (e.g., rational points) to reconstruct the flow. 

Although a trajectory representation for a monotone flow T with profile 
Z is not unique, there is a special representation U(x,to,ti) that is right- 
continuous in x G [0, Z(to)] for every t\ > to: 

^ x ' t0 ' tl) = lx, x>z(t ). (12) 

The concrete way of defining U(x,to,ti) for x > Z(to) is inessential for 
our purposes, and we often will simply ignore points (x,to,ti) with x > 
Z(io). 

It is often convenient to understand a monotone flow as a triple T = 
(r, Z, U), where Z is the profile of T, and U is one of the trajectory repre- 
sentations of r. 

We are now ready to define a solution of our main equation ([9]). 
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6 Solution of the SPDE 



In this section, we define a solution of ([9]), prove its existence, uniqueness, 
and the Markov property. 

Let (f2, J-, P) be a complete probability space with a filtration (J-t)t>o sat- 
isfying the usual conditions (right-continuity and completeness), and let W 
be a Brownian sheet on w.r.t. (J r t)t>o, i-e., W is a centered Gaussian ran- 
dom field indexed by Borel subsets of , such that for any Borel set A C 
with finite Lebesgue measure \A\, (W([0, t] xA)) t >o, is an (J r t) t >o-martingale 
with (W([0, •] x A))t = t\A\. We refer to [Wal86] for more background on 
martingale measures. 

We define a solution to © on [0, T] as a random monotone flow r : f2 — > 
M[0,T] such that with probability 1, the triplet (T,Z,U) has the following 
properties for some trajectory representation U of T: 

1. For any to > and any ii > to, T* '* 1 is measurable w.r.t. Tt x - 

2. For any x > and (to-,t\) G A, U(x,to,ti) is measurable w.r.t. J^i- 

3. For any to > 0, any x > 0, almost every cj € {x < Z(to)}, and all 
*i > to, 

U(x,t ,t 1 )=x+ — ' ' dt+ / / l[o,^(* I t ,t)](!/)W'(dyxdt). 

(13) 

Let us recall that the idea of this definition is to define a natural existence 
and uniqueness class of solutions by means of consistent co-evolution of 
monotone graphs and trajectories of individual points. 

Theorem 5 For any filtered probability space (Q, J-, (J-t)t>o, P) satisfying 
the usual conditions and a Brownian sheet W w.r.t. (J-t), there is a unique 
solution (r*°'* 1 )(t 0) t :L ) e A of equation ([9]) in the above sense. 

Proof: Let us construct a solution first. We begin with Z. Let us start 
with an auxiliary equation 

Z(x,t) = x + t+ [ [ l[ ,z(x, s )](y)W(dy x ds) 

JO JR 

for x > 0. Mimicking the SDE methodology, it is easy to construct a unique 
strong solution of this equation defined up to a stopping time 

T n = inf 1 1 > : Z(x, *)<-!. 
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Letting n — > oo and noticing that is a no-exit singular point for the dif- 
fusion Z{x, •), we can extend this solution to the whole time semi- axis R+. 
Using the uniqueness, we see that with probability 1, for all m € N and 
t > 0, 

Z(l/m,t) > Z(l/{m + l),t), 

so that the limit Z(t) = lim m _ i , 00 Z(l/m, t) is well-defined. It is easy to 
check now that Z{t) is a strong solution of equation 

Z(t) = t+ f [ l [0 , zi s)](y)W(dy x dt). 

JO JR 

It is also easy to check that this strong solution is a unique one by mimicking 
the proof of the uniqueness theorem of Yamada and Watanabe for SDEs, 
see Proposition 2.13 in [KS881 Chapter 5]) 

For any two nonnegative rational numbers x, t, we define a stochastic 
process (U(x, to, ti)) tl >t as follows. If x > Z{to), we set C7(x,to,ti) = x; if 
x < Z(to), we set U(x, to, •) to be the unique strong solution of equation (fl~3j) 
on [to, oo) (the uniqueness can be established easily using methods borrowed 
from the theory of SDEs). 

For each t± > to the resulting map x i-> U{x,to,t\) is nondecreasing 
on [0, Z(to)]: if x\ < X2 then U(xi,to,ti) < U(x2,ta,ti). In fact, if the 
opposite inequality holds for some t\ > to then U(xi,to,t') = U (x2,to, t') for 
some t' € [£o 5 il]j an d the uniqueness implies that U{x\,to, ■) and U(x2,to, •) 
coincide after t' . 

Now, for any (to,t\) E A with rational to, the (random) map x — >■ 
U{x,to,t\) defined above uniquely determines a monotone graph r io,il . Us- 
ing the continuity of trajectories U(x, to, t{) and their monotonicity property, 
it is easy to show that r* '* 1 depends continuously on rationals to and t\. 
Therefore, we can extend T* '* 1 by continuity to all values {to,t\) G A. 

Knowing monotone functions U{x,to,t\) for rational values of x and to, 
and all t\ > to allows to define 

U(x, to, ti) = mi{U(y, t, h) : t, y € Q, t < t , y € [0, Z{t)], U(y, t, t ) > x}, 

for all other values of to and x £ [0, Z(to)], and it is easy to verify that thus 
defined random field U satisfies the requirement of the Theorem. 

To prove the uniqueness, it suffices to notice that the result of each step 
in the above construction of is a.s.-uniquely defined. □ 

The proof of the theorem provides a construction of a unique solution of 
SPDE ©. In the next section, we study some properties of the solution. 
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Remark 3 The set of all pairs of rationals x,t in the proof of Theorem^ 
can be replaced by any countable and dense set in R^ 

The following statement follows directly from the Markov property of 
U(x,t ,t{) in tr- 

Lemma 6 (Markov property) For any to > 0, (T to,tl ,ti > to) is a Markov 
process. 

7 Discontinuities of the solution 

The reason for introducing monotone flows the way we did in SectionsEJandE] 
is the presence of discontinuities of solutions with respect to x € [0,Z(t)\. 
If not for these discontinuities our analysis would have been easier, and in 
this section we will show that they are, in fact, intrinsic to the solution, thus 
justifying our choice to work with monotone graphs. 

Let us call a straight line segment connecting two points (x\,yi) and (#2,2/2) 
on the plain vertical if x\ = x%. 

Lemma 7 Let < to < t\. Then, with probability 1, r* '* 1 contains nonde- 
generate vertical segments, i.e., any monotone function representing r* '* 1 
is discontinuous, with shocks associated to these vertical segments. 

Proof: To prove this lemma we need the following standard result (see e.g. 
the proof of Theorem 26 in |Pro041 Chapter 2] for this result in the context 
of realizations of stochastic processes) . 

Lemma 8 If f : [0, z] — > R is a bounded variation function then its quadratic 
variation Q(f) defined by 

n-l 

Q(f) = lim J2 (/(*(* + l )l n ) ~ f( zi / n )f 

n— ¥00 — * 

i=0 

satisfies: 

Q(f)= E (/(*+) -/(*-)) 2 > 

xeA(f) 

where and A(/) is the set of all discontinuity points of f. In particular, f 
is continuous on [0,z] iff Q(f) = 0. 
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Let us compute the quadratic variation of U(-, to, t±). For n € N we 
introduce (omitting dependence on n) 

k 

x k = -Z(t ), k = 0,...,n, 
n 

and 

Vk(ti) = U(x k ,t ,ti) - U(x k -i,t ,ti), k = l,...,n. 
Ito's formula implies 

dV k 2 (h) =2V k {h)dV k {h) + V k (h)dh 

= J°u i r d *i + 2Vfc(*i)W([l7(xfc_i, to, ti), t , ti)] x dti) 
+ Vfc(*i)rfti. 

Let 

Qn(tl) = ^l 2 (tl) + ... + K 2 (il)- 

Then Q n (t„) = (Z(t )) 2 /™, and 

„ ((l) .« +2 r'%w rfi 

" A Z(t) 

+ 2 2Z / V k (t)W([U(x k ^,t ,t),U(x k ,t ,t)]xdt) 
k=i Jt ° 

+ / Z(t)dt. 

Jt 

Let us define Q(t) = lim^oo Q n (t) and v = inf{t > to : > 0}. If 

v > to, then taking the limit as n — > oo in both sides of the equation above 
at ti = v, we see that all terms converge to zero except for Z{t)dt. To 
obtain the convergence to zero for the martingale stochastic integral term, 
it is sufficient to see that the quadratic variation in time given by 

In.. V ,. v n 

(zZ ^(W[^fc-i,io,t),C/(x fe ,t ,t)] xdt)) = / ^2v£(t)dt, 
\k=i Jto I v Jt ° k=i 

converges to zero as n — > oo. 

Since Z(t)dt is a strictly positive random variable that does not de- 
pend on n, we obtain a contradiction which shows that v = t , so that 
Q(t) > for any t > t . □ 



16 



8 Convergence 



This section is the central part of the paper. Here we prove that the discrete 
infinite random genealogy tree r converges in distribution in an appropriate 
sense under appropriate rescaling to the continuum random tree given by 
the monotone flow described in the previous section. 

8.1 Trees as monotone flows. Main result. 

To start with, we introduce a monotone flow associated with a realization 
of the infinite discrete tree r. This procedure is analogous to that of linear 
interpolation for discrete time random walks leading to Donsker's invariance 
principle. 

We begin with an imbedding of the tree in the plane. Recall that there 
are X n > 1 vertices in the ra-th generation of the tree. For i £ {1, . . . , X n }, 
the i-th vertex of ra-th generation is represented by the point (n, i — 1) 
on the plane. The parent-child relation between two vertices of the tree 
is represented by a straight line segment connecting the representations of 
these vertices. 

Besides these "regular" segments, we shall need some auxiliary segments 
that are not an intrinsic part of the tree but will be used in representing the 
discrete tree as a continuous flow. Suppose a vertex i in n-th generation has 
no children. Let j be the maximal vertex in generation ra + 1 among those 
having their parents preceding i in generation n. Then an auxiliary segment 
of type I connects the points (ra, i — 1) and (ra + l,j — 1). If vertex 1 in ra-th 
generation has no children, points (ra,0) and (ra + 1,0) are connected by an 
auxiliary segment of type II. Auxiliary segments of type III connect points 
(n, i — 1) and (ra, i) for 1 < i < X n — 1. 

Every bounded connected component of the complement to the union of 
the above segments on the plane is either a parallelogram with two vertical 
sides of length 1, or a triangle with one vertical side of length 1. One can 
treat both shapes as trapezoids (with one of the parallel sides having zero 
length in the case of triangle). 

For each trapezoid, we shall establish a bijection with the unit square 
and define the monotone flow to act along the images of the "horizontal" 
segments of the square. A graphic illustration of the construction is given 
on Figure HJ and we proceed to describe it precisely. 

Each trapezoid L of this family has vertices go,o = { n , ^o,o)> 50,1 = ( n > ^0,1)1 
5i,o = (™+Mi,o),5i,i = (n+l,ii,i), where z ,i-«o,o G {0, 1} and n,i-ii,o G 
{0,1}. 
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Figure 1: Construction of the continuous monotone flow. 



Then, for every a G (0, 1) we define 

9m( a ) = 9m,o + a(g m ,i ~ 9m,o), m = 0,1, 

and 

g L (a, s) = g^{a) + s(g[(a) - g^(a)), s G (0, 1). 

This definition introduces a coordinate system in L, i.e., a bijection be- 
tween L and the unit square (0, l) 2 . We are going to use it to construct 
the monotone map associated with the tree for times to,ti assuming that 
there is n G {0} U N such that n < to < t\ < n + 1. Let us take any x 
such that (x, to) belongs to one of the trapezoids L. Then there is a unique 
number a(x,to) G (0,1) such that g L {a{x, to), {to}) = x, where {•} denotes 
the fractional part. We can define g t °' tl (x) = g L (a(x , t ) , {h}) . This strictly 
increasing function can be consistently and uniquely extended by continuity 
to points x such that (x,to) belongs either to a regular segment in the tree 
representation or an auxiliary segment of type I or II. This function g* '* 1 
also uniquely defines a monotone graph r* '* 1 = r*°'* 1 (r) depending contin- 
uously on to, ii- Next, if we allow to and t\ to take values n and n + 1, then 
we can construct the associated monotone graph as the limit in (M, p) of 
the monotone graphs associated to the increasing functions defined above 
(as to — > n or t\ — > n + 1). Notice that the resulting monotone graphs may 
have intervals of constancy and shocks (i.e., contain horizontal and vertical 
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segments). Now we can take any (io^i) € Aoo and define 

fio,*i = f [ti],ti Q f[h]-l,ti Q _ _ _ Q p[t ]+l,[to]+2 fto,M+l_ 

which results in a continuous monotone flow (r* '* 1 (T))(t ,t 1 )eA oo • 

To state our main result we need to introduce a rescaling of this family. 
For every n G N, we define 

r^ il (r) = {f — ,-^) r^^Gf^Hr)), (io^GA^. (14) 
[\pn pnj J 

For each T > we can consider the uniform distance pt on monotone flows 
in M[0, T] introduced in (jlip . and define the locally uniform (LU) metric 
on M[0, oo) via 

oo 

d(r 1 ,r 2 ) = ^2~ m ( / ) m (r 1 ,r 2 )Ai). 

m=l 

Theorem 6 As n — > oo, the random field (Tn ' ll (t)) ^ t(htl ^ e ^ converges in 
distribution in LU metric to the stochastic flow of monotone graphs solving 
equation ([9]). 

The rest of the section is devoted to the proof of this theorem. First, it is 
sufficient to prove the convergence in distribution in p m for all m. We take 
m = 1 without loss of generality. Due to the classical Prokhorov theorem it 
is sufficient to demonstrate two facts: 

1. The sequence of distributions of r n (r) is tight w.r.t. p\. 

2. Any limit point for the sequence of distributions coincides with the 
distribution of the flow generated by SPDE. 

Subsection 18.21 is devoted to the proof of the first statement and subsec- 
tion 18.31 to the second one. 



8.2 Tightness 

We have to show that for any e > there is a set K C M[0, 1] that is 
compact in d\ and P{r n G K} > 1 — e for all n (from now on we do not 
distinguish between T n and its restriction on [0, 1].) 

The construction below depends on the values of constants a > 0, j3 > 
7 > and b > 0. For m G N, consider a set 

R m = {(&2-"\ j2~ m ) : < k < 2 m ; < j < 2 m ( 1+Q )} , 
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and also define R = {J m R m - For a function / we denote its /3-H61der 
constant by Hp(f): 

^ m |/(*i)-/(* 2 )| 

tl,t 2 :ti^2 Fl — *2| P 

Let if m be the set of all monotone flows (T, Z, £/") in M[0, 1] such that 
sup t Z(t) < 2 ma , Hp(Z) < 62" 17 , and for each (t ,x) G i? m satisfying 
x < Z(to), Hft(U(x,to,-)) < 62™t. Denote K = f)™= mo K m . Note that 
we suppress the dependence of K on b and mo- 

Lemma 9 For any 6 and m , the set K is precompact in M[0, 1]. 

Proof: For any sequence (F n , Z n ,U n ) ne ?$ in .KT, we must show that it 
contains a convergent subsequence. 

Using the Holder continuity, the Arcela-Ascoli compactness criterion, 
and the classical diagonal procedure, we can extract a subsequence n' such 
that for each m G N and each each (t ,x) G R m , U n /(x,to,-) converges 
uniformly to a limiting trajectory Uoo(x,tQ, •) with 'H/3(U 00 (x,to, •)) < 62 m7 , 
and Z n i converges uniformly to a limiting trajectory Z^. 

We need to show that there is a unique continuous monotone flow T, 
compatible with these trajectories and that T n > converges uniformly to V 
Let us construct = {Too, Z^,, U^). Notice that Z^ is already at our 
disposal. Take any point (to, x) G [0, 1] x R + such that x < Z^to). For any 
m, we set 

k = k(m) = [2 m t ], 
j-=j-(m) = [2 m (x-2 m i-(2- m f)], 



oo 
oo ■ 



and 

Then, 

and 



j+ = j+(m) = [2 m (x + 2™ 7 ■ (2- m f)} + 1. 
j+ 2 -m _ j- 2 ~™ < 2(2~ m + 2™ 7 • 2~ m/3 ), 



■- 2 -ra + 2 m 7 . <X< j+2" m - 2" 17 • (2~ m f\ 

For sufficiently large m, j+2~ m < Z 00 (k2~ m ), so that U oo (j ± 2- m ,k2- m ,t ) 
are already defined. Using the bound on the Holder constant, we get 

x - 2 - m - 2 ■ 2 m b-ft < U oo {j-2-' m ,k2~ m ,t ) < x 

< U oo {j + 2~ m ,k2~ m ,t ) < x + 2~ m + 2 ■ 2 m (~<-P\ 
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For any t\ > to, we define 

U-(x,t ,t 1 ) = S upU oo (k(m)2- m ,j-(m)2- m ,t 1 ), 

m 

and 

U + (x,t ,t 1 ) = m{U oo (k(m)2- m ,j + (m)2- m ,t 1 ). 

m 

It is clear that U~(x, to, ti) < U + (x, to, ti) and if x < y then U + (x, to, h) < 
U~(y,to,ti). In particular, both quantities are monotone in first argument 
and there is a monotone graph T^,'* 1 such that U~(x, to, ti) and U + (x, to, ti) 
are its left continuous and, respectively, right continuous representatives. 
Also any monotone function compatible with (U(x, to, ti), (to, x) G R), must 
represent the same monotone graph. 

The resulting family of monotone graphs (T^'* 1 ) is easily seen to satisfy 
all the properties in the definition of a monotone flow. Let us prove only 
Property 2, the continuity of T^'* 1 in to and t\. It is sufficient to show that 
T^,'* 1 converges to the graph of the identity map on [0, (t)] as to and ti 
converge to t from below and above respectively. To see the latter, we imbed 
a small time interval (to,ti) into some dyadic interval [k2~ m ,(k + l)2~ m ] 
with a large m. On {k2~ m , (k + l)2~ m ], the monotone flow displaces all 
dyadic points with denominator 2~ m by at most 2 m ^~^ . Since 2 m ^~^ — )■ 
as m — > oo, the continuity follows. 

Now we shall prove the convergence to this monotone flow T^. For any 
e > 0, we will show that for sufficiently large values of n', c?i(r oo , r n /) < e. 
First, for any m G N, we can find no(m) such that for all n' > rio(m) and 
all (to, ^) € Rmi \\U n (x,to, ■) — U^x, to, Oil < e /2, where || • || denotes the 
sup-norm. For any (to, a;) choose points k,j~,j + as above and denote 

^ = ^(^2—^2—,^), 
yt = U O o(j ± 2- m ,k2- m ,t ). 



Then 

x - 2~ m - 2 ■ 2 m{ ^ <y-,<x< y+ <x + 2~ m + 2- 2 m{ ^ , 
x - 2~ m - 2 ■ 2 m ^~^ <y^<x<y+<x + 2~ m + 2- 2 m ^~^. 

The image of x under T^ ,'* 1 (viewed as a multivalued map) for every t\ > to 
is contained in the image of [y^nUn^ under T^ ,'* 1 . Therefore, the definition 
of y^t and y^ implies that if n' > no(m), then any point in the image of 
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[y~,, y„,] under T^ /* 1 is at most at e/2 from the image of [y~ , y+] under T^'* 1 . 
On the plane, the distance between any point (x, x') G T^ ,'* 1 and the graph 
r^'* 1 restricted to does not exceed 2~ m - 2 • 2 m ^- /3 ) + e/2. So 

it is sufficient to choose m so that 2~ m — 2 • 2 m ( 7_/3 ) < e/2 to have that 
distance bounded by e. Similarly, the same bound holds true for the distance 
between any point (x,x') G T^,'* 1 and the restriction of T^ / 1 to [y^iVoo]- 
We conclude that for the above choice of m and for any n! > no(m), the 
uniform distance between T n > and Too is bounded by e, and the convergence 
implying the precompactness of K follows. □ 

Our next goal is to show that for any e > 0, the numbers a, j3, 7, b can 
be chosen so that P{r n G K} > 1 — e for all n. 

Lemma 10 Suppose there are positive numbers C, r such that for all m £ N, 
all (to,x) G R m , a continuous process U(x,to,-) satisfies 

E\U(x,t ,t) -U(x,t ,t')\ r < C\t-t'\ r / 2 , t <t<t'. (15) 

If 

r/2-l>r/3, (16) 
then there is a constant C±(C,f3,r) such that for all m G N and any b > 0, 

P{ sup Hf,(U(x,t ,-))>bY> m \<C 1 (C,P,r) . 

[(t ,x)eRm ) 

Proof: Let us fix x,to and (for brevity) denote U(t) = U(x,to,t), 

Let us estimate P{T-L/3(U) > c} for a c > 2. If %p(JJ) > c, then there are 
times t < t' such that 

\U{t')-U(t)\ >c{t'-tf. 
We can find m such that 2-' m < t' - t < 2- m+1 . Then 

\U(t')-U(t)\ >c2" m/3 . (17) 



There is an integer j such that j2 m G [t, t']. We can find sequences («&)£ 
and (n' k ) k x L m such that each and /{' fc is or 1 and 

00 

k=m 



00 

m 



22 



j2-m - t = ^K fc 2- fc . 
k=m 

Assuming that on each respective dyadic interval of length 2~ k , the incre- 
ment of U does not exceed (1 — 2~^)c2~^ k ~ l in absolute value, we have 



117(0 " U(t)\ < (1 - 2-/ 3 )c2~ 1 £ 42"^ + £ K fc 2" 



-8k 



< (1 - 2~ /3 )c2" 1 • 2 ^ 2 

< c2- m P, 



k=m k=m 
oo 

j3k 



which contradicts (fT7|) . Therefore, denoting z/(/3) = (1 — 2~ /3 )2~ 1 , we see 
that for some k there is a dyadic interval [I2~ k , (/ + l)2 _fc ] of length 2~ k such 
that 

\U{l2~ k ) - U((l + l)2- fc )| > v{l3)c2-P k . 



We use Markov's inequality and condition (115]) to see that an upper bound 
for the probability of this event is: 

oo 

^2 fc p{|[/(Z2- fe ) -U((l + l)2- fc )| > v{p)dT pk \ 



k=0 

_. oo 

— — 2 k E\U(l2- k ) - U((l + l)2- k )\ r 2^ k 
' [P> k=0 

C O k -kr/2 rl3k < Ci(C,f3,r) 



k=0 



where 



C ' 



k=0 

Notice that C\{C,j3,r) < oo due to (fT6j) . The lemma follows as we set 
c = 62 7m since there are 2^ 2+a ^ m points in R m . □ 

The proof of the following lemma is similar to that of Lemma [10\ and 
we omit it. 

Lemma 11 Ifr/2 — 1 > r/3 and Z is a continuous process satisfying 

E\Z(t') -z{t)\ r < c\t'-t\ r / 2 , M'e[0,l], (18) 
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then 

ty—r-ym 

?{Up{Z) > b2^ m } < C^C^r)—^. 

Lemma 12 Let 

2 + a<rj. (19) 

Suppose r = (r, Z, U) is a monotone flow such that processes U satisfy the 
conditions of Lemma \l(A and process Z satisfies the conditions of Lemma\ll\ 
Suppose there is a constant C such that 

E sup Z(t) < C. (20) 
te[o,i] 

then there are b = b(a, ft, 7, r, C, C) and mo = wio(a, ft, 7, r, C, C) such that 

P{r g k} > 1 -£. 

Proof: 



P{reif}<^Pi sup Hf)(U(x,to,-)) >b2^ m 

m=l 



(to,x)eR„ 
00 



m=l 



+ P i sup Z(t) > 2 am \ . (21) 

Lemma [10] and condition (|19p imply that for sufficiently large 6, the first 
term in the r.h.s. is less than e/3. It also follows from Lemma [11] that for 
sufficiently large b, the second term in the r.h.s. is less than e/3. Finally, we 
can use (|20p and Markov's inequality to choose mo so that the last term in 
the r.h.s. is less than e/3 completing the proof. □ 

Since we want (fT6|) and (fT9|) to be satisfied along with 7 < ft, we need 
to have 

2 + a<rft < 

This can be satisfied if we take r = 8, 1/4 < 7 < /3 < 3/8, < a < 87 - 2, 
and our goal is to show that all other conditions are satisfied with this choice 
of parameters for the sequence of random monotone flows T n = (T n , Z n ,U n ) 
defined in (fl"4"|) . 

From now on we assume that U n is the special right-continuous trajectory 
representation of T n , see (fl~2j) . 
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Lemma 13 For all n, the process Z n restricted to times 0,n~ , 2n~ , 3n~ , . . 
is a submartingale with respect to its natural filtration. For any r £ N, 
EZ^(t) is uniformly bounded in t £ [0, 1] and n £ N. 

Proof: It is sufficient to consider the process X n counting vertices in the n- 
th generation of the discrete random tree and show that it is a submartingale 
and that for every r there is a constant c r such that for all n£N, 

EX r n < c r n r . (22) 

The latter follows from the following one: for every r there is a polynomial 
P r —i of degree r — 1 such that 

E[X r n+1 \X n ]<X r n + P r ^(X n ). (23) 

Taking expectations of both sides of (f23j) and applying a straightforward 
induction procedure in r and n, we obtain the growth estimate (I22p , 

To prove (|23|) . we define -P(r, /c) to be the set of all partitions of r into k 
nonnegative integers: 

P(r, k) = {(ji, ...,3k)- Ji > J2 > • • • > ifc > 0, ji+j2 + ...+ jk = r}, 



■■ii k Ph---Pi k 
■i{ k Ph---Pi k 

(24) 

Here ii, . . . , denote the number of children of each of k vertices of gener- 
ation n, and in the third line we used the symmetry of the expression and 
grouped together monomials producing the same partition j € P(r + l,k) 
under the monotone rearrangement of their degrees (we agree that 0° = 1). 
The numbers B m have been introduced in ([2]). The positive constants Kj(r) 
satisfy 

£ Kj(r) = k r +\ (25) 

jSP(r+l,fc) 
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and write 

E[X r n+l \X n = k] = l (ii + --- + ik) r+1 Pii 



k 

0<h,...,ik<D 



r + 1 



n l 2 



k o<h,...,i k <D j Vii,--->ifc 

0<ii,...,i fe <DjgP(r+l,fc) 

i E Kj(r)B jl B h ...B jk . 

jeP(r+l,k) 



Notice that if k > r + 1 then 
A" 



:.i 1.0.0...0 = k(k - 1) ... (A; - r) = - P r (fc), 

where P r (k) is a polynomial of degree at most r. In this case, due to ([25]) . the 
sum of all other constants equals P r (k). Now (|23p follows straight from (|24p . 
The fact that (X n ) is a submartingale follows from 

E[^ n +i — X n \a(Xo, . . . ,X n )} = u, 

which was established in [Bak]. This identity can also be derived from the 
computation in the proof of Lemma [TU □ 

Lemma 14 The conditions of Lemma [T3 o- re satisfied for U n uniformly 
in n, if we choose r = 8, l/4<7</3< 3/8, < a < 87 — 2. 

Proof: We have to prove that (fT5|) is satisfied for U n (x,t , •) uniformly in 
n,x,to,t,t'. Suppose first that (to,x) is a grid point, i.e., t^n and xfrn are 
integers. Denote by Vj the size of progeny in generation j generated by first 
xfj^n + 1 vertices in ton-th generation. Then U n (x,to, j /n) = (Vj — l)/(/m). 
Let us compute 



V, 



m+l 



1 

link 
1 

link 
1 



Vrn 

an 



I X m 
un' un 



k 

fin 



^2 (ii + ■ ■ ■ + k)(h + • • • + ik)Ph ■ ■ -P 
Q<h,...,i k <D 

h 

[IB 2 + l(k - 1)] 



hi2Ph-..pi k 



Therefore, 



ank 
/in V k 



m+l 
un 



Vrn 

/in 



Vrn 

fin 



J_ x m 

lin un 



k_ 

un 



ll_ 

n k 



(26) 



Since < l/k < 1, the process A n [t) interpolating between the values of 

[nt]-l 

A n {t) 



[nt\-l 



Yl 

un — ' Xn 

p 3=0 



3 
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is a 1-Lipschitz process, and 

B n {t) = U n (x,t ,t) - A n (t) 

restricted to t = mn^ 1 , m = 0,1,..., n, is a martingale. Since A is 1- 
Lipschitz, we now need to prove that B n satisfies a moment estimate anal- 
ogous to (|15p with constant in the r.h.s. independent of n. If t = mn -1 and 
t' = m'n~ l , then due to Burkholder's inequality (see ShifMi. Section VII. 3]) 
and Minkowski's inequality (sec [Shi96, II.6] ) , there is an absolute constant 
C 2 such that 

4 



'm'-l 



E\B n (t') - B n (t)\ s < C 2 E \J2 i B n((j + - BnijrT 1 )) 

\j=m 

ml — 1 

£ {BM + l^-Bnijn- 1 )) 2 



< C 2 



m'-l 



<c 2 [Yl IK^ca+i)"" 1 )-^"- 1 )) 2 



j=m 



< C 2 (m' - m) 4 max (B n ((j + l)n _1 ) - ^(j^ 1 ))' 



< C 2 n\tl - i) 4 maxE |S n ((j + I)™" 1 ) - £„(jn 

To estimate the r.h.s., we use the following statement: 
Lemma 15 There is an absolute constant C3 such that 

E |5 n ((j + l)n- V £n(jO)| 8 < % 
The desired moment estimate follows directly: 

E\B n (t') - B n (t)\ 8 <C 2 C 3 (t> -t) 4 . 
Proof of Lemma [T5l 



-1M 8 



B r , 



i + i 



n 



n 



1 

< — + 
/m 



V, 



1"; 



/in /in 



so that, by convexity, 

'i + i 



z?„ ( - 

n 



< 2' 



(jm) 1 



+ 2 7 E 



/in /in 



and the lemma follows directly from the next result: 
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Lemma 16 For any even r > 2, there is a constant c' r such that for all n 
and all j, 

^.m-^rl. (27) 

\ n n J n J 
Proof: Let us estimate E(Vj+i - Vj) r . 

E[(v j+1 -v j y\v j = i,x j = k] 

=- ^2 (n + ••• + «/ - l) r (h + ... + ik- k)p h ...p ik 

ii,...,ifc 

+ E (h + --- + ii- l) r Ph ■■■Pi k 

h,—,ik 

= 1 E -!) + ■■■ + (*»- l)) r ((»i - 1) + . . . + (»* - 1)K ---Pik 

h,...,ik 

+ E ((«l-l) + --- + (»i-l)) r Pii---Pi* 



i ^ ((<i - 1) + . . . + (i, - l))^ . . . 



ll,...,lfc 



+ r E ((*i-i) + --- + (*i-i)) r ((ii+i-i) + --- + (*fc-i)K--. 



+ E ((«i-i) + --- + («i-i)) r Pii---Pi* 

ii,...,ifc 



We begin with 53. Since 



Pii+i---Pi k = 1 i 

k+i, — ,*fc 



^3= E E (n - ■■■(^ 

j€P(r,l)h,-,ii 

= E K jBh---B jv 

jeP(r,l) 

where all monomials with similar degree j are grouped together, and 
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Notice that B\ = 0, so that if one of the terms in j\ + . . . + ji = r equals 1, 
the contribution from that term to 53 is 0. In other words, only partitions 
of r containing no l's contribute to S3. Each of these partitions contains 
at most r/2 nonzero components, and therefore the associated coefficient 
Kj (the number of monomials associated to this partition) is bounded by 
a degree r/2 polynomial in The same analysis shows that S2 = 0, and 
S\ is bounded by a polynomial of degree r/2. We conclude that there is a 
constant C' r (independent of k) such that 

E[(^'+i " V 3 ) r \ V, = l,Xj =k]< C' r rl\ (28) 

The estimate (j22j) implies now that 

m+i - VjY < c'rtxf < J T fi\ 

so that ((27} holds true. □ 

We assume now that the desired moment estimate for {to,x) on the 
grid holds true with a constant C. We need is to estimate E\U n (x,to,t') — 
U n (x,to,t)\ 8 for arbitrary n, x, to, t, t' . 

Let j = [tn], f = [t'n] and N = j' — j. We shall assume that ton < j; 
the proof can be easily modified to treat the opposite situation. 

We are going to consider three cases: (i) N = 0, (ii) N = 1, (iii) N > 2. 

In case (i), the evolution of U n (x,to,t) is linear in t S \jn , (j + l)n~ l ], 
so that, comparing it to the evolution along the regular edges of the tree 
imbedding, we see that 

f U n (x,t ,t')-U n (x,t ,t) \ 8 _ ( Un(.x,t ,^)-U n (x,to,^) \ 8 
V t'-t ) l/n ) 

<n 8 £ E [(U n (x, t ,^)- Un {x, t , m Unix, *g] , 

1 

+ E[(V j+ i - V 3 f\V 3 = l]P{u n (x,t , i) € (£, ^]}) • 
Inequality ([28} implies now that we can continue the estimate as 

E(U n { X ,to,t') -U n {x,to,t)f 

< 2(f - tfC> £(i + 1) 4 P {Unix, to, i) € (jt, %]} 

I 

< 2(t' — t) 8 /i 4 n 4 Cg E (U n (x, t , + ^) 4 

< 2(1/ - t)YC 8 Ct, (29) 
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where 



CI = sup sup E I Z n (t) H 

neNte[o,i] V V n 



P 4 



and in the last line we used the fact that t' — t < n . In case (ii), 

'u n (x,t ,t') - u n (x,t ,ty 



< E 



t'-t 

'U n (x,t ,2±k)-U n (x,t ,t)\ fu n {x,t ,t')-U n (x,t ,^) 



< F ( Un(x,t ,^)-U n (x,t ,l) \ 8 / Un ( x ,t ,^)-U n (x,t ,^) \ S 

and for each of the terms in the r.h.s. we can proceed as in case (i) and 
obtain an estimate analogous to ([29]) . 

Case (hi). Let us estimate E(U n (x,to,t') — U n (x,to,t))+ first. Let £ = 
\U(x,to,t)fin\. Suppose U(x,to,t') > U(x,to,t). Then the monotonicity of 
the flow implies that at least one of the following two conditions holds true: 

L U fa, i, i±i) - U(x, t , t) > \{U{x, t , t') - U(x, t ,t)). 
2. Forr)=[U(x,to,i£)»n] + l, 

Consequently, 

(U(x, t ,t') - U{x, t ,t)) + < 2{R 1 AR 2 A R 3 ), 



where 



Ri 


= u\ 


' '_L 






Ri 


= u\ 


/ f] 






\fin 


Rz 


= u\ 


f 7] 









U(x,t ,t), 



it l fin n 



Therefore, 



///; // // J fin n 



E(U(x,t Q ,t') - xf + < 2 8 (ERf + ERl + ER%) 
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Notice that the random variables R\,R2, R3 are defined in terms of the 
flow trajectories emitted from lattice points. Although these lattice points 
defined through £ and rj are random, they are measurable with respect to 
the filtration of the flow, and therefore the above estimates imply that there 
is a constant C such that 

E(U(x, t ,t') - U(x, t , *))+ < C ^ ~ *) 4 - 

We can estimate E(U(x,to,t') — U(x,to,t))^_ similarly, which completes the 
proof. 

8.3 Characterization of limit points 

We begin with two auxiliary lemmas. 

Lemma 17 Let mo € N. For each n € N, take a random number V(n) of 
vertices in generation mon measurable with respect to the history of the tree 

(n) 

up to generation mon. For m > mon we denote by V m the total size of the 
progeny of these vertices in generation m. (In particular, Vmon = V(n).) If 

P{V(n)/n > n~ 7 } ->■ 0, n -> 00, 

for some constant 7 S (0, 1), independent of n, then, for any m\ > m , 

V {n) p 

Vm P „ 

sup > (J, n — > 00. 

mon<m<min ^ 

Proof: Let X m ,m > is the process of total population sizes. If we show 
that 

V {n) p 

sup — > 0, n — > 00, (3U) 

mon<m<min -A-m 

then the lemma will follow since we can write 

T/O) T/0 1 ) V 

Vm Vm -^m 
sup = sup — — • sup , 

mon<m<min W mon<m<min -A-m m(,n<m<m\n fl 

and use the convergence in distribution of sup m()ri < m < miri — — to the maxi- 
mum of a diffusion process. 

To prove (f30j) . we define a sequence of events 

A n = {X mon > n 1 -^ 2 ; V(n) < n 1 ^}, 
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and for any e > write 



Yn 



V„. 



P sup -^ >e \<P{A c n ) + P[l sup -^ >£ \nA n 

Clearly, P(-A^) — >• as n — > oo due to the assumptions of the lemma and 
the fact that the distribution of X n /n converges weakly to a distribution 
with no atom at 0. 

'y(n) ' 
* m 



For the second term, we notice that y-j^j is a nonnegative bounded 

supermartingale w.r.t. filtration (J-m), where for each m = itiqu, . . . ,min, 
(n) 

Tm is the sigma-algebra generated by Xi, Vi, I = mow, . . . , m. In fact, 



V. 



(n) 

+1| y(n) — TV — U 
L m+1 

1 h + . . . + ij 

k ^— ' i\ + . . . + ik 

0<ii,...,ik<D 
ii + ...+j fc ^0 



(«! + ... + i k )p h ...p 



\ £ (*! + ••• + • • • Pi k - rPo 1 



0<h,...,i k <D 



< 



I 1 

k ~ k 
I 

k' 



Po 



Doob's maximal inequality for nonnegative supermartingales implies 

({ V \ \ n 1-7 In 1 '"'/ 2 
P I \ sup >£ nij < - >0, n -> oo, 

and (I30p is proven. □ 

Lemma 18 Zei 7 G (0,1/2). For any x > 0, t > 0, there is a constant C 
such that for all n£N, 

p{\U n (x,t, Mtl) - «| > n" 7 } < Cn-^ 

Proof: It is easily seen that for large n, due to the monotonicity of the 
flow, if 

[nt] + 1 \ [xfin] 



U n I x , t , 



> n 
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then 



and, moreover, 



v ■ ([xri [nt] + l\ _M > 
\ fin n ) an 



fin 



Also, for large values of re, if 



L x/m] [ntj + l^ [x/ire] 



/in 



> re 



-7 



7? 



£/ n I x,t, 



nt] + 1 \ [x/jn] 



then 



and, moreover, 



^ / [xfin] f [nt] 



< — n 

fin 

+ 1 \ [x/jn] 



-7 



< 



/ire 



Therefore, 



/[x/m] [nt] + l^ [x/m] 
u n I , t. 



n 



fin 



> re~ 7 - -. 
n 



< P{\U n ( 

The moment estimate (|28p from Lemma [TBI implies that 



> n 



"7 _ 2 



n \ /in ' ' n 



[x/m] 



/in 



,H ,x 
fi z n z n 



and the desired estimate follows from Markov's inequality. □ 

Now we proceed to prove that any limiting point for the sequence of 
distributions of monotone flows T n has to coincide with the monotone flow 
solving the SPDE Q as discussed in Section [6j 

Consider (xo,i,io) = (0, 0), a sequence of times < t\ < t<i < . . . < tj., 
and for each i = 1, . . . , k, a sequence of nonnegative numbers 



xa< ... < x U (i)- 
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For each i = 0, . . . , k and every j = 1, . . . , l(i) we define a family of processes 

[nU] 



WW 



-131 



{xjjfmj [ntj] + l 



un 



t < 
'), t> 



[ntj 



Due to Lemmas [T7] and [TBI it is sufficient to show that as n — > 00, 
these processes jointly converge in distribution in sup-norm to a nonnegative 
diffusion process (U^) with drift 

11' A 

b lJ (y,t) = -^-l {t>u . o< yi j<y 01 }, i = 0,...,k, j = 1,... ,Z(i), 
and diffusion matrix 



= (y hjl ^yi2j2)l{t>t il VU 2 ; 0<y tljl ,y t2J2 <yoi}' 



h,h = 0, ...,k, j 1 = l,...,l(h), j 2 = l,...,l(i 2 ). 

For any i and any y, these coefficients are constant on t £ [ti, ti+i), they are 
continuous and bounded in y, and define a well-posed martingale problem. 
Let 

K J (t) = U^([nt]/n) 
Ift>t,and^(M)<C/„ 01 (M), 



So, we define 



[nt] + 1 I I F 
n n ' 



Tit] 



w n 



M\ 1 



(31) 



and 



M%(t) = U*>(t)-B%(t--). 

n 



Next, a simple calculation based on the martingale property of M„ and (|3T1) 
shows that if f > U, Ut jl ff) < U^ff), and Ut h ff) < U* 1 ^) then 



n n ' 



hh ( N + 1 




M + 1 



n 



nt]/n 



r;.^|^ + -)i.F 



1 



5 2 -l 



nt]/n 



(32) 



II l Mil /' I:"/. 



[nt] 



1 ) I 1 h • 



1 S2-1 



j? 



n jjm 



nt] 
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Processes Un J1 and Un J2 describe the evolution of sizes of two subpopu- 
lations, i.e., they are rescaled versions of vertex-counting discrete processes 
and V( 2 \ To compute the first term in the r.h.s. of (|32p . we take 
< h < fa < k and write 

y(2) | y (l) _ , y (2) _ , y 
V m+l V m+l\ V m ~ l li v m ~ l 2,^m — « 

: T ^ (*1 + ■ ■ ■ + %)(*1 + • • • + ih)( i i + • • • + i k)Pi 1 ■ ■ -Pi k 
ii,...,i k 

= i(ii5 3 + - 1) + h{fa - 1) + /i(/ 2 - 1))B 2 + /i(/ 2 - l)(fc - 2)). 



Denoting 



- ,,„,£/« ( M 
n 



and making a simple calculation, we can rewrite the r.h.s. of (|32H as 

/1 fl + R(k,fa] 



n \n kn 
where R(-, •) is a uniformly bounded function. We denote 
[ni] 

^n Jll2J2 (*) = - X/ 1 / TO >nft i .,Vfe„1:C/: ii HM^ £ 7:2i2 f MKC/i?ifMA X ( 33 ) 



n 

m- 



^ A {m>n(t il Vt ia );Ui lJl (^),I^ , * ! (^)<^ 1 ( J ^)} 

x([/^(^) A ^ 2 (^))x (34) 



x 1 + 



n^ 01 (f) 



The above calculation shows that 

^WK l31 (t)-^ J112,2 W 

is a martingale. 
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We are going to use Theorem 4. 1 from |EK86l Chapter 7] . Identities (f3Tj) 
and (|33p imply two main conditions of that theorem: 



sup 

0<t<l 



b^(U(s),s)ds 



and 



sup 

0<t<l 



(U(s),s)ds 



0. 



The other conditions of Theorem 4.1 }EK86] are concerned with the jumps 
of processes A and B, they can be easily proven using the moment estimates 
above. The theorem allows us to conclude that the process U n converges 
weakly in the Skorokhod topology on D[0, 1] to the continuous diffusion 
process with drift and diffusion given by b and a respectively. In fact, The- 
orem 4.1 is given in |EK86j for time- independent coefficients, and here we 
invoke its straightforward time-dependent generalization. 

Convergence of U n in distribution in Skorokhod topology to a contin- 
uous process implies convergence in distribution in sup norm to the same 
process. Moreover we can conclude that the U n , the linear interpolated ver- 
sion of U n also converges in distribution in sup norm, which concludes the 
demonstration of the theorem. 



9 Connection to superprocesses 

There is an important connection of our results to the theory of super- 
processes. Superprocesses are measure-valued stochastic processes describ- 
ing the evolution of populations of branching and migrating particles, see 
e.g. |Daw93]. The limiting SPDE that we have constructed is similar to 
the genealogy in the Dawson-Watanabe superprocess with no motion con- 
ditioned on nonextinction, see [Eva93j . [EP90] . |DK99j . 

Our approach is more geometric then the superprocess point of view. For 
the superprocess corresponding to our situation, the continual mass momen- 
tarily organizes itself into a finite random number of atoms of positive mass 
(corresponding to discontinuities of the monotone maps in our approach). 
The mass of these atoms evolves in time analogously to equations ([7]), but 
our approach helps to understand what happens inside the atoms by unfold- 
ing the details of the genealogy giving rise to a contentful cocycle property 
whereas in superprocesses the cocycle action essentially reduces to evolving 
the masses of atoms. The monotone flow we construct can be initiated with 
zero mass and just one particle whereas in superprocesses one has to start 
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Figure 2: Stochastic foliation constructed for 600 generations of a tree. 



with positive mass. The related historical superprocesses do describe the ge- 
nealogy, but they require nontrivial particle motion and ignore the ordering 
and monotonicity issues that are central to our approach. 

Notice that the continuous monotone flow we construct, describes a fo- 
liation of the random set {(t, x) : < x < Z(t)} into diffusion trajectories. 
The geometry of this stochastic foliation is suprisingly complicated. Due to 
the presence of shocks in monotone maps (corresponding to the atoms of 
superprocesses), the stochastic foliation cannot be obtained, say, as a con- 
tinuous image of the foliation of a rectangle in "horizontal" segments, and 
the flow is very far from a flow of diffeomorphisms. 

We hope that our results are interesting from the point of view of graph 
theory since they describe what a typical large tree looks like. In fact, 
combining the results of this paper with those of Section [2j we conclude that 
a typical embedding of a large ordered rooted tree in the plane if rescaled 
appropriately looks like a stochastic foliation described by SPDE (fT3|) . It 
would be interesting to obtain rigorously a direct convergence result that 
would not involve the intermediate infinite discrete tree. However, currently 
this kind of result is not available. 

Figure [2] shows a realization of a pre-limit monotone flow for a large 
random tree. Every tenth generation is split into about ten subpopulations, 
their progenies are tracked and shown on the figure. 

References 

[AP98] David Aldous and Jim Pitman. Tree-valued Markov chains derived 
from Galton- Watson processes. Ann. Inst. H. Poincare Probab. 
Statist, 34(5):637-686, 1998. 



37 



[Bak] Yuri Bakhtin. Thermodynamic limit for large random trees. 

Accepted at Random Structures Algorithms. Also available at 
http://arxiv.org/abs/0809.2974. 

[BH08] Yuri Bakhtin and Christine Heitsch. Large deviations for random 
trees. J. Stat. Phys., 132(3):551-560, 2008. 

[BH09] Yuri Bakhtin and Christine Heitsch. Large deviations for random 
trees and the branching of RNA secondary structures. Bull Math 
Biol, 71(1):84-106, 2009. 

[Daw93] Donald A. Dawson. Measure-valued Markov processes. In Ecole 
d'Ete de Probabilites de Saint-Flour XXI — 1991, volume 1541 of 
Lecture Notes in Math., pages 1-260. Springer, Berlin, 1993. 

[DK99] Peter Donnelly and Thomas G. Kurtz. Genealogical processes for 
Fleming- Viot models with selection and recombination. Ann. Appl. 
Probab., 9(4):1091-1148, 1999. 

[EK86] Stewart N. Ethier and Thomas G. Kurtz. Markov processes. Char- 
acterization and convergence. Wiley Series in Probability and 
Mathematical Statistics: Probability and Mathematical Statistics. 
John Wiley & Sons Inc., New York, 1986. 

[EP90] Steven N. Evans and Edwin Perkins. Measure-valued Markov 
branching processes conditioned on nonextinction. Israel J. Math., 
71(3):329-337, 1990. 

[Eva93] Steven N. Evans. Two representations of a conditioned superpro- 
cess. Proc. Roy. Soc. Edinburgh Sect. A, 123(5):959-971, 1993. 

[Ken75] Douglas P. Kennedy. The Galton- Watson process conditioned on 
the total progeny. J. Appl. Probability, 12(4):800-806, 1975. 

[Kes86] Harry Kesten. Subdiffusive behavior of random walk on a random 
cluster. Ann. Inst. H. Poincare Probab. Statist., 22(4):425-487, 
1986. 

[KS88] Ioannis Karatzas and Steven E. Shreve. Brownian motion and 
stochastic calculus, volume 113 of Graduate Texts in Mathematics. 
Springer- Verlag, New York, 1988. 

[LG99] Jean-Frangois Le Gall. Spatial branching processes, random snakes 
and partial differential equations. Lectures in Mathematics ETH 
Zurich. Birkhauser Verlag, Basel, 1999. 



38 



[Pro04] Philip E. Protter. Stochastic integration and differential equations, 
volume 21 of Applications of Mathematics (New York). Springer- 
Verlag, Berlin, second edition, 2004. Stochastic Modelling and 
Applied Probability. 

[Shi96] A. N. Shiryaev. Probability, volume 95 of Graduate Texts in Math- 
ematics. Springer- Verlag, New York, second edition, 1996. Trans- 
lated from the first (1980) Russian edition by R. P. Boas. 

[Wal86] John B. Walsh. An introduction to stochastic partial differential 
equations. In Ecole d'ete de probabilites de Saint-Flour, XIV— 
1984, volume 1180 of Lecture Notes in Math., pages 265-439. 
Springer, Berlin, 1986. 



39 



